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^ \ Abstract 

We searched for redshifted O emission lines from the possible warm-hot intergalactic medium (WHIM) 
f~>) ■ surrounding the cluster of galaxies A2218 at z = 0.1756 using the XIS instrument on Suzaku. This cluster 

\^ \ is thought to have an elongated structure along the line of sight based on previous studies. We studied 

systematic uncertainties in the spectrum of the Galactic emission and in the soft X-ray response of the 
' detectors due to the contamination building up on the XIS filters. We detected no significant redshifted O 

lines, and set a tight constraint on the intensity with upper limits for the surface brightness of O VII and 
O VIII lines of 1.1 x 10~ 7 and 3.0 x 10 -7 photons cm -2 s _1 arcmin~ 2 , respectively These upper limits are 
significantly lower than the previously reported fluxes from the WHIM around other clusters of galaxies. 
> ■ We also discuss the prospect for the detection of the WHIM lines with Suzaku XIS in the future. 
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■ 1. Introduction 

_ _ i 

Based on several N-body simulations of cosmological large-scale structure formation (e.g., Cen & Ostriker, 1999; 
Dave et al., 2001; Chen et al., 2003), a significant (30-50%) fraction of baryons is thought to reside in the form of 
gas in a 'warm-hot' phase (T = 10 5-7 K), which is hard to detect with the instruments currently in operation. This 
warm- hot gas, whose density is 10 -6 -10~ 4 cm~ 3 , is called the warm- hot intergalactic medium (WHIM). The firm 
detection of the WHIM is important because it is the most promising candidate for the "missing baryons"; i.e., it is 
thought to explain the discrepancy between the baryon density observed in the local universe (Fukugita et al., 1998) 
and that in the distant universe (Rauch et al., 1997) or that calculated from the observed fluctuations of the cosmic 
microwave background (Spergel et al., 2003). 

The WHIM may be detected via emission or absorption lines from highly ionized elements in UV or X-ray spectra. 
Given the relative elemental abundances, — on the order of 0.1 to 0.3 of Solar — oxygen is the most promising clement 
to provide detectable transitions at the expected temperatures. Over 40 systems show Ovi absorption features, based 
on observations from FUSE and HST (e.g. Danforth & Shull, 2005). However, their analysis shows that the baryon 
density in the temperature range T = 10 5 ~ 6 K is about an order of magnitude lower than the expected WHIM level. 
Therefore, most of the missing baryons seem to reside in a hotter phase at T = 10 6 " 7 K, which can be probed with 
X-rays. There is only one sightline with a likely detection of the WHIM via OVII or Ovill absorption lines so far 
(Nicastro et al., 2005). Their results suggest that most of the missing baryons can be explained by the absorbing 
clouds inferred from their data. However, recent work by Rasmusscn et al. (2006) and Kaastra et al. (2006) showed 
that the detection with the Chandra LETGS by Nicastro et al. (2005) was not consistent with RGS observations of 
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the same sight line. Thus, the existence of the WHIM in the hotter phase is not yet confirmed. 

Several possible detections of WHIM emission have been reported. Kaastra et al. (2003) and Kaastra (2004) report 
the detection of redshifted O VII and O VIII emission lines as well as soft X-ray (E < 0.5 kcV) excesses in the spectra of 
7 clusters of galaxies out of 21 measured systems. Finoguenov et al. (2003) also reported detection of O VII and OVIII 
emission lines in the outskirts of the Coma cluster. Fujimoto et al. (2004) and Takei et al. (2006) report not only 
emission lines in the spectrum of cluster outskirts but also absorption lines in the spectrum of a quasar behind the 
clusters, although the significance of the absorption measurement is not high (< 3d). Fujimoto et al. (2004) observed 
OVIII features in the Virgo cluster, while Takei et al. (2006) detected NelX in the Coma cluster. They estimated 
the path length of absorbing/emitting medium assuming that the absorber and the emitter are the same cloud, and 
concluded that the length scale of the medium exceeds the virial scale of the clusters. There are counterarguments 
against the discovery of O VII or O VIII emission lines associated with clusters. For example, Bonamentc et al. (2005) 
and Bregman & Lloyd-Davies (2006) concluded that some of the emission lines observed by Kaastra et al. (2003) are 
likely due to field-to-field variation of the soft X-ray background or Galactic emission. Such controversies occurred 
because the CCD instruments cannot clearly distinguish O lines emitted by the gas associated with clusters of moderate 
redshift (z < 0.05) from those of Galactic origin. These arguments suggest the importance of careful target selection 
and accurate measurement of background level. Note that the observations of redshifted emission lines also give 
important clues to solve another mystery "cluster soft excess", the excess emission in E < 0.3 kcV, reported since the 
1990s in the extreme ultraviolet and the X-ray bands from EUVE and ROSAT (e.g., Lieu et al., 1996a; Lieu et al., 
1996b; Mittaz et al., 1998; Bowyer et al., 2004). Thus, the study of X-ray emission in cluster outskirts brings us 
important information about the cluster system itself and its environment. 

Since there are only a few reports of possible emission lines from the WHIM, further sensitive study of other systems 
is desirable. In this paper, we report our observation of A2218 with the XIS (Koyama et al. 2006) onboard Suzaku 
(Mitsuda et al. 2006). A2218 is a well-known cluster because of its strong gravitational lensing arcs. Its redshift is 
0.1756 (Struble & Rood, 1999), where one arcminute corresponds to 179 kpc with the cosmology we have adopted. 
The cluster is thought to be still in a dynamically young state — X-ray analysis suggests that the cluster core is 
not in hydrostatic equilibrium due to an ongoing or recent merger (Pratt et al., 2005; Machacek et al., 2002) in the 
line-of-sight direction. The distribution of galaxy velocities based on the optical study of Girardi et al. (1997) indicates 
substructure in the line of sight, and Smith et al. (2005) studied the mass structure in detail and have concluded that 
A2218 is unrelaxed. The dynamics of the cluster is considered to cause the discrepancy between the mass estimated 
from gravitational lensing and that determined from X-ray observations; gravitational lensing indicates ~ 3 times 
higher mass than X-ray estimation for the central (~ 20") region. These observations indicate an elongated structure 
in the line-of-sight direction and suggests existence of a large-scale filament in this direction. If this is the case, we 
would expect a high surface brightness of the WHIM emission. The fairly large redshift of A2218 would also help us 
observe a significant energy shift of emission lines with the XIS instrument. Thus, A2218 seems to be one of the most 
suitable targets to search for the WHIM emission from Suzaku. 

We assume a Hubble constant of 70 km s _1 Mpc -1 or h?o = 1 and fi m = 0.3, f^A = 0.7. The solar metal abundances 
used in this paper are given by Anders & Grevesse (1989). Unless otherwise stated, errors in the figures are at the 
68 % confidence level and at the 90 % confidence level in the text and tables. 

2. Observations 

We carried out four observations to study the warm-hot gas in the A2218 vicinity with the XIS instrument onboard 
Suzaku: two on the cluster and two in offset regions whose locations are different from each other. The offset 
observations were performed to measure the foreground Galactic contribution, which can produce a large ambiguity 
in the study of soft emission around clusters. Fig. 1 is a map of the ROSAT R4 band after removing bright point 
sources (Snowden et al., 1997). This energy band is sensitive to OVII and OVIII emission lines, which are strong in 
the Galactic emission, and hence is suitable to study their fluctuations. 

We examined the ROSAT map (Fig. 1) and selected the location of offset pointings so that the Galactic emission 
level is similar to that in the vicinity of A2218. Also, the distance from A2218 is large (> 1° or 10.7 Mpc at the 
A2218 redshift), so the data from the offset regions should be free from the putative large-scale filament around 
A2218. The observed locations are shown in Fig. 1 and Table 1, as well as the position of A2218. Hereafter we call 
the four observations A2218-1, A2218-2, Offset-A, and Offset-B, respectively. A2218-1 and Offset-A were observed in 
early October, 2005, and the others in late October. The observation sequence numbers and dates are summarized in 
Table 2. 

The XIS is an X-ray CCD camera, which consists of three front-illuminated (FI) sensors and one back-illuminated 
(BI) sensor. The camera covers the 0.2-12 keV energy band with an energy resolution of 130 cV at 5.9 kcV and 
40 eV at 0.5 keV. Its field of view is a square of 18' x 18'. The BI and FI sensors have different advantages: the 
former has higher quantum efficiency in the soft energy band, while the latter shows a lower internal background 
level. An important characteristic of the XIS is the absence of a large low-energy tail in the pulse-height distribution 
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function, even in the very soft energy (E < 0.5 keV) band. Those characteristics give the XIS an excellent capability 
to study the low-energy emission lines, such as from the WHIM. After the launch of Suzaku, a time and position 
dependent contamination of the XIS optical blocking filter (OBF) was found. The source is probably outgassing from 
the satellite. The level of contamination increases with time and is different from sensor to sensor. The time and 
position dependence of the contamination thickness has been empirically modeled by the XIS team. At the time of 
our observations, the effective area at 0.5 keV was on average 25% lower than that of the pre-launch calibration. See 
Koyama et al. (2006) for a fuller description of the instrument. The two observations of A2218 were 25 days apart, 
and the count rates at 0.5 keV were different by 9% and 13% for the BI and the combined FI sensors, respectively. 
This is consistent with the expected drop of 8% for both sensors within statistical errors. 

The XIS instrument was operated in normal mode in all the observations presented in this paper. We used event 
files of version 0.7 products 1 . Events of 3 x 3 and 5x5 observation modes were combined. The first few kiloseconds 
of each observation, when the pointing direction had not yet stabilized, were excluded. The resulting image of A2218 
in the 1.0-4.0 keV energy band is shown in Fig. 2. Data from the two observations and the four CCDs are combined. 
The white circles indicate the region where we extracted the spectra. The radii are 3' and 8' from the cluster center, 
corresponding to 540 and 1430 kpc at the source, respectively 
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Fig. 1. ROSAT R4 band image around A2218 (Snowden et al., 1997) in Galactic coordinates. The field of view of XIS in the four 
observations are indicated by black squares. Most of the emission in this band is O VII and O VIII line. 



Table 1. Position of the A2218 and offset pointing 



(RA, Dec) (I, b) 

A2218 (16 h 35 m 54 s , 66°13'00") (97.7449, 38.1235) 

Offset-A (16 h 17 m 48 s , 65°27'36") (97.7423, 40.1239) 

Offset-B (16 h 39 m 31 s , 66°13'31") (96.4000,38.1002) 



Version processing is an internal processing applied to the Suzaku data obtained during the SWG phase. Aspect correction and fine 
tuning of the event time tagging are skipped in this version. Times when the elevations from the bright and dark Earth are > 20° and 
> 5° are excised in standard data processing. 
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Table 2. Suzaku observations of A2218 and offset pointing 





Sequence 


Date 


Exposure 


Net exposure 


Net exposure 




number 




after processing 


using all COR 


using COR> 8 Gev c" 1 


A2218-1 


100030010 


2005-10-01-2005-10-02 


46.4 ks 


44.9 ks 


38.2 ks 


A2218-2 


800019010 


2005-10-26-2005-10-27 


32.8 ks 


32.3 ks 


28.8 ks 


Offset-A 


100030020 


2005-10-02-2005-10-03 


44.6 ks 


44.6 ks 


39.0 ks 


Offset-B 


800020010 


2005-10-27-2005-10-27 


15.0 ks 


15.0 ks 


12.0 ks 




Fig. 2. Image of A2218 in the 1.0—4.0 keV band observed with the XIS. The events of the two observations obtained by four CCDs 
were summed. The spectra were extracted from the annulus between the two white circles. Vignetting effects were not taken into 
account for this image and background events were not subtracted. The grid indicates the coordinates in J2000.0. 

3. Spectral analysis 

3.1. Analysis method 

We removed flickering CCD pixels, which cause a large noise component below 0.5 keV. We also excluded apparent 
point sources that were found in the Offset-A image. Although the fluxes of these point sources were lower than 
the detection limit for the A2218 image, we removed them because we arc interested in the upper limits of the soft 
excess emission. The spectra were then extracted in the annular region between 3' and 8' from the cluster center 
using XSELECT distributed in HEASOFT 6.0.4. We define the cluster center as (DETX, DETY) = (499, 530) in 
detector coordinates for all observations and all sensors. Note that the boresight of all XIS sensors coincide within 20". 
We have chosen the inner radius to be 3' in order to exclude the bright central region, in which the strong emission 
of the hot ICM hampers the study of warm-hot emission. We also excluded the region outside of 8' because the 
position dependence of the contamination on the OBF is not well known there. Although the XIS sensors are always 
illuminated by 55 Fe calibration sources at the field edge, we did not remove these regions in the present analysis. 

To estimate and subtract the internal background, we also extracted spectra from night Earth observations dis- 
tributed by the XIS team. It is known that the variation of the internal background spectrum is well correlated with 
the cut-off rigidity (COR) at the position of the satellite; the smaller the COR value is, the larger the background level 
becomes. To estimate accurately the internal background level, we collected events when the detector was looking 
at the dark Earth (night Earth events) and sorted them by COR. We extracted the spectrum for each 1 GeV cT 1 
interval of COR, and then added them weighted by exposure time for the respective COR range in the actual obser- 
vation. This process gives different night Earth spectra for the four observations, since the detailed distribution of 
COR was different among the observations. The background spectra were extracted from the same region in detector 
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Fig. 3. Spectra of A2218 taken with the (left) BI and (right) FI chips, where three FI spectra were averaged. Black, red, green 
and blue points show A2218, Offset-A, Offset-B, night earth spectra, respectively. The strong peaks at 5.9 keV and 6.5 keV in the 
right panel are due to the calibration sources. 

coordinates as the corresponding observation in order to avoid possible systematic effects due to positional variation 
of the detector background. 

We used response matrix files (RMFs) ae_xi[0123]_20060213(c).rmf distributed by the XIS team for spectral fitting. 
On the other hand, we generated the ancillary response files (ARFs) with the arf builder 'xissimarfgen', which is based 
on ray-tracing (Ishisaki et al. 2006). The change of the effective area with time in the soft X-ray energy range, due 
to the increase of OBF contamination, was taken into account by the following method. The ARFs were created 
separately for each observation, according to the observation date. We used the ae_xi[0123]_contami_20060525.fits 
contamination tables to model the composition and position dependence of the contaminant, in which the C/O ratio 
of contaminant was assumed to be 6. Note that the effective area for diffuse sources depends on the sky distribution 
of the flux, because the quantum efficiency varies with the detector position due to vignetting and contamination 
thickness. We adopted the Suzaku XIS image (1.0-4.0 keV) of A2218 for the incident flux distribution, in which the 
two observations by the four sensors were all summed (see Fig. 2). We assumed a uniform flux distribution for the 
offset pointings. 

In the analysis shown below, the spectra and ARFs for A2218-1 and A2218-2 were added using the ftools mathpha 
and addarf, respectively. Further, the spectra and response (RMFs and ARFs) of the three FI sensors were combined 
using ftools mathpha and marfrmf, respectively. We confirmed that consistent results were obtained with this treatment 
within the statistical error. We fitted spectra of the offset observations in the 0.35 keV < E < 5.0 keV energy range for 
the BI and 0.40 keV < E < 5.0 keV except the energy of anomalous response at the Si K-edgc, 1.825-1.840 keV, for 
the FI. When we fitted the spectra of A2218, the energy range was extended up to 7.0 keV in order to cover the Fe-K 
lines. We excluded 5.85-5.95 keV and 6.45-6.55 keV, because the slight difference in intensity of Mn K emission lines 
from the 55 Fe calibration sources between A2218 observations and night Earth observations causes relatively large 
residuals. The spectra of A2218 (black), Offset-A (red), and Offset-B (green), as well as the night Earth spectra (blue) 
are shown in Fig. 3. Left and right panels are for BI and FI spectra, respectively. The difference in the detector area 
due to point-source exclusion was corrected in the figure for illustrative purposes. We mainly performed simultaneous 
fits for the two spectra obtained by the BI and the combined FI sensors. However, it was also confirmed that the 
spectrum with the BI sensor only gave consistent results. This indicates that the uncertainty in the FI response around 
0.5 keV is not crucial compared with the statistical errors in the BI spectrum. 

Evaluation of the systematic uncertainties is crucial in constraining the emission from the WHIM. Firstly, the 
emission is in the soft X-ray region where the XIS has various systematic uncertainties, so we have to look into many 
possible effects. Secondly, the observation was carried out only three months after the launch of Suzaku, and time 
variation of the detector response needs to be carefully considered. 

The first uncertainty is the detector response. The largest systematic uncertainty in the detector response is the 
thickness of contaminant on the OBF of the XIS. In the 0.4-0.6 keV band, most of the photons we observe are thermal 
emission from the ICM. If we underestimate the contaminant thickness, the flux from the ICM would be overestimated 
in E ^5 1 keV, because the temperature and abundance are strongly constrained by the data in the higher energy range 
E>1 keV. This naturally leads to an underestimation of the soft-excess flux. Therefore, we also considered a thicker 
contaminant model at the upper limit of the uncertainty range. We generated ARFs with 20% thicker contaminant 
using xissimarfgen assuming the observation occurred seventeen days later than it actually did. This 20% increase, 
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which exceeds the uncertainty quoted by the XIS team (~ 10%; Koyama et al. 2006), is a reasonable value, given the 
additional uncertainties in the response to diffuse emission and in early observations. The energy resolution of the 
XIS worsens with time. At the time of our observations, the energy resolution had degraded by ~ 20 eV (FWHM) 
at 5.9 keV (from 130 eV to 150 eV). This effect was investigated by smoothing the response with a Gaussian. Since 
we have no information about the degradation at ~ 0.5 keV, we tried Gaussians with sigmas of 5 eV, 10 eV, 15 eV, 
20 eV, 30 eV and 35 eV. With this smoothing, the spectral resolution at ^0.5 keV increases from 40 eV to 42 eV, 
46 eV, 53 eV, 62 eV, 81 eV and 91 eV, respectively. The results are presented later. The internal background of 
the XIS is quite stable, and there were no "background flares" during our observations. However, when COR is low 
(COR< 6 GeV c _1 ), the internal background level becomes comparable to the X-ray background in the E > 3 keV 
range. Therefore, we also extracted spectra with the condition COR > 8 GeV c _1 to examine this systematic effect. 

Another uncertainty is the spatial variation of the Galactic emission. The ROSATKA map shows ~ 10% variations 
among the A2218, Offset- A and Offset-B fields. The purpose of the two offset observations (Offset- A and Offset-B) 
was to look at the spectral variation in the Galactic emission and to include it in the analysis. Besides the inclusion 
of this variation, we also extended our upper limit for the soft excess by assuming 10% fainter Galactic emission. 

3.2. Offset pointings 

Before the analysis of the A2218 spectrum, we analyzed data from the offset pointings to estimate the spectrum of 
the Galactic emission. The diffuse X-ray background at high Galactic latitude can be divided into three components: 
i.e. the local hot bubble (LHB), the Milky Way halo (MWH) and the extragalactic power-law (CXB) components 
(Snowden et al., 1998). Typical temperatures are about 0.1 keV for the LHB and 0.2-0.3 keV for the MWH, while the 
CXB spectrum has a photon index of 1.4. The LHB is thought to surround the solar system with a ~ 100 pc scale, 
and hence it has no Galactic absorption. The CXB component is extragalactic and known to be uniform over the 
entire sky. In contrast, the level of LHB and MWH components vary from position to position. Lumb et al. (2002) 
reported that the mean deviation of the 0.2-1.0 keV intensity is ~ 35% from field to field. 

Our purpose here is to constrain the spectra of the two Galactic components so as to obtain reliable background data 
for the estimation of the warm emission around A2218. Spectra for an annular region with r = 3' — 8' were produced 
for the Offset- A and Offset-B data. XSPEC vll.3.2r was used in the spectral analysis. The spectra, from which non 
X-ray background was subtracted, were fitted with a model for the sum of LHB, MWH and CXB. CXB and MWH 
components were absorbed by Galactic absorption characterized by Ah = 3.24 x 10 20 cm -2 (Dickey & Lockman, 1990). 
The collisionally ionized thermal plasma model APEC was used in XSPEC to fit the LHB and MWH spectra, with 
abundance and redshift fixed to 1 solar and 0, respectively. The photon index of the power-law CXB model and the 
temperature of the LHB were fixed to 1.4 and 0.08 keV, respectively. The normalizations of the FI and BI sensors 
were allowed to take different values. 

The best-fit parameters and the model curves are shown in Table 3 and Fig. 4. In Table 3, parameters except 
for the FI/BI ratio are for the BI spectra. Acceptable fits were obtained for the two spectra and the parameters of 
the Galactic emission were well constrained. The normalization of the CXB component agrees with previous reports 
(Kushino et al., 2002). The flux in the 0.5-2.0 keV band of Offset-A and Offset-B observation differs by ~ 10%. Not 
only the flux, but also the shape of the best-fit model is different between the two observations, which is either due to 
the poor statistics of the Offset-B observation or real variations in the Galactic emission. 



Table 3. Best-fit parameters of offset pointings 



Parameter Offsct-A Offset-B 

LHB kT (keV) 0.08 (fixed) 0.08 (fixed) 

LHB Normalization 3.2 ± 1.0 x 10~ 6 6.2±^ x 10~ 6 

MWH kT (keV) 0.158t£;5j 2 ! 0.248^;^ 

MWH Normalization 21 6.7±J'jj x 10~ 7 5-0±\i x 10" 7 

CXB Photon index 1.4 (fixed) 1.4 (fixed) 

CXB Normalization 13 8.5 ± 0.5 x 10~ 7 8.5 ± 0.7 x 10~ 7 

FI/BI ratio 0.84 ±0.05 1.04 ±0.09 

X z /dof 154.03/149 100.45/93 

~ J n c n H oIV/4:1t(Da{1 + z))' z per solid angle in units of 
10 14 cm -5 arcmin -2 , where n e is the electron density, nn the 
hydrogen density, and £*a the angular size distance. 
b In units of photons cm" 2 i 
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Fig. 4. Spectra and best-fit models, and residuals of offset observations, (left) Offset-A, (right) Offset-B. Black and red represents 
BI and FI sensors. 

3.3. A 22 18 

3.3.1. Single temperature model 

The A2218 spectrum for the radial region of 3' — 8' was analyzed here. The summed FI and the BI spectra were 
fitted simultaneously. Before testing the warm emission in the model, we first fitted with a model only including the 
thin thermal plasma (ICM) and the background emission determined from the offset-region spectra. We tried two 
models for the background: namely, the best-fit models for the BI spectra determined from the Offset-A and Offset-B 
observations. We fixed all the parameters for the background spectrum consisting of LHB, MWH and CXB at the 
best-fit values. The redshift of the hot ICM was fixed at 0.1756. We allowed the FI/BI normalization to be free 
again. The results of the fit are summarized in Table 4 and in Fig. 5. The different background (Offsct-A or B) gives 
an ICM temperature differing by about 10%, while similar metal abundance values were obtained in the two cases. 
The temperature and abundance are generally consistent with previously determined values; XMM-Newton spectra 
indicates kT ~ 5 keV in 3' — 5' and Z = 0.13 ± 0.04 solar in the central 5'.1 region (Pratt et al., 2005), while Chandra 
estimated kT to be 6.9 ±0.5 and Z to be 0.20 ±0.13 in the central 5'.1 region (Machacek et al., 2002). The detailed 
analysis of the temperature and abundance of the ICM will be reported elsewhere. 

The spectra were well fitted with this simple model, suggesting that no obvious soft excess is present. The residual 
of the spectral fit should indicate a feature, such as that caused by redshifted oxygen lines, if warm emission indeed 
gives a significant contribution. Both panels in Fig. 5 suggest some features around 0.57 keV and 0.68 keV: the former 
is a positive residual, while the latter negative. Both features are near the strong Galactic O emission lines: O VII at 
574 eV and O VIII at 654 eV, and not identical to those of redshifted O lines at the redshift of A2218. Their levels 
are less than that of Galactic emission by a factor of 3-10. Hence, they may originate in the fluctuation or incomplete 
modeling of Galactic emission. Note that the O VIII line (654 eV) at the source is redshifted to 556 eV, which is close 
to the strong OVII line (574 eV) in the Galactic emission. Thus it is rather difficult to detect the warm gas with this 
line. The redshifted OVII line (574 eV shifts to 488 eV) is also close to the peak caused by the oxygen edge of the 
detector deadlayer and OBF, in particular for the FI sensors (see Fig. 4 of Koyama et al. 2006). These situations 
make the XIS sensitivity to warm gas for a source with redshift similar to that of A2218 poorer than it is at other 
redshifts. We will carry out a quantitative evaluation in the next section. 

3.3.2. Constraint on redshifted O lines 

We performed another spectral fit by adding two Gaussian emission lines to the model, in order to test for the 
existence of redshifted O lines. The energies of the lines were fixed to 488.22 eV and 555.99 eV, which correspond 
to O VII (resonance) and OVIII lines at z = 0.1756. Although the energy of the OVII line could be at most 10 eV 
lower due to the contribution of intercombination and forbidden lines, wc confirmed that the results shown below are 
not affected by that difference in energy. The intrinsic widths were fixed to zero. The temperature and abundance of 
the ICM component were free parameters, while all the background parameters (for LHB, MWH and CXB) including 
the normalizations were fixed at the best-fit values determined with the Offset data described in § 3.3.1. The free 
parameters were the temperature kT, abundance Z and normalization of the ICM, the surface brightness I of redshifted 
OVII and OVIII lines, and FI/BI normalization ratio. We again tried the two background models, determined from 
the Offset-A and Offsct-B observations, respectively. 

The best- fit parameters and improvement in x 2 values are shown in Table 5. The obtained surface brightness / 
of the lines is small and consistent with zero; the improvement in \ 2 1S l ess than 2.71 (90% significance for 1 free 
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Table 4. Best-fit parameters of A2218 with single temperature model 



Parameter 



Offset-A Case 



Offsct-B Case 



LHB kT (kcV) 
LHB Normalization a 
MWH kT (keV) 
MWH Normalization 3 
CXB Photon index 
CXB Normalization 13 



0.08 (fixed) 
3.2 x 10 -6 (fixed) 

0.158 (fixed) 
6.7 x 10 -7 (fixed) 

1.4 (fixed) 
8.5 x 10 -7 (fixed) 



0.08 (fixed) 
6.2 x 10 -6 (fixed) 

0.248 (fixed) 
5.0 x 10 -7 (fixed) 

1.4 (fixed) 
8.5 x 10 -7 (fixed) 



ICM kT (keV) 


5 40+ u - 2V 


6.00 ±0.22 


ICM Z (solar) 


0.20 ±0.04 


0.21 ±0.04 


ICM z 


0.1756 (fixed) 


0.1756 (fixed) 


ICM Normalization 3 


2.08 ±0.03 x 10 -5 


2.00 ±0.03 x 10 -5 


FI/BI ratio 


0.89 ±0.01 


0.90 ±0.01 


X'fdof 


475.21/475 


520.72/475 



j n c nn cIV/4it(Da(1 + z))' z per solid angle in units of 
10 14 cm -5 arcmin -2 , where n c is the electron density, nn the 
hydrogen density, and Da the angular size distance. 



In units of photons cm 



keV arcmin at 1 keV 



[Vol. , 




1 2 
Energy (keV) 




Energy (keV) 



Fig. 5. Spectra fitted with the best-fit models, and residuals of A22f8 observations, in which the data from the two observations 
were merged. Galactic components are fixed to (left) Offset-A, and (right) Offsct-B models, respectively. Black and red indicate BI 
and Fl sensors. The contribution of the Galactic emission, determined by offset observations, are also indicated by dashed lines. 



parameter). We then constrained the allowed O line intensity range. The upper limits for the O VII and O VIII / were 
estimated from that causing an increment of the x 2 value by 4.0 (2a limit) over the minimum (best-fit) level. In this 
process, we allowed the other free parameters (i.e., the ICM temperature, abundance, normalization, and FI/BI ratio) 
to vary. The upper limits of the two lines were determined separately, and the values are also shown in Table 5 in 
parentheses. Fig. 6 shows the model with the upper-limit intensities of the O VII and O VIII lines in the upper 2 panels; 
the left panel is the case with the background of Offset-A, while the right panel is for Offset-B. The upper limits are 
shown with solid lines, while the background spectrum is shown with dashed lines. The upper limits of the two O 
lines are 10-20% of the background level at the center energy of each line. 

It should be noted that the O VIII line is also produced by the ICM (T ~ 5 keV), while O VII is not. The intensity, 
which depends on the temperature and O abundance, is 2 x 10 -8 photons cm -2 s -1 arcmin -2 assuming the parameters 
in Table 4. However, this value contains an uncertainty; while the abundance was mostly determined by the Fe K 
and L features in the fit, Fe and O abundances may have a different value. We estimated the upper limit of I for 
the redshiftcd OVIII line from the WHIM in the extreme case of no O in the ICM, because this situation gives the 
largest contribution from the WHIM. The upper limit then increased to 1.6 x 10 -7 photons enr 2 s -1 arcmin -2 and 
1.8 x 10 -7 photons cm -2 s -1 arcmin -2 for the Offset-A and Offset-B backgrounds, respectively. 

Next, we investigated how the upper limits of I change by considering the systematic uncertainties described in § 3.1. 
We already checked one of these uncertainties in Table 5: variation of the Galactic emission between Offset-A and 
Offset-B fields. The difference is < 20%. When we used the ARFs with 20% thicker contaminant, the upper limits of the 
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VII and O VIII lines increased by ~ 130% and ~ 30%, respectively. These limits further increase by ~ 15% and ~ 30% 
for the respective lines, if we adopt 10% fainter Galactic emission. On the other hand, the upper limits do not increase 
when we select the data for the low background condition (COR> 8 GeV c _1 ) or when the energy resolution of the 
detector was artificially degraded by 5-35 eV. To summarize, the conservative upper limits are derived by assuming no 
O in the ICM and employing the ARFs with 20% thicker contaminant and 10% fainter Galactic model for the Offset-B 
observation; the values are 1.1 x 10 -7 photons cm -2 s _1 arcmin -2 and 3.0 x 10 -7 photons cm -2 s -1 arcmin -2 for 
OVII and OVIII lines, respectively. We adopt these values as the upper limits of the O line intensity around A2218, 
considering both statistical and systematic errors. The spectra with the model with these upper limits are shown in 
the lower panel of Fig. 6. Note that even after considering these systematic errors, the upper limits are still ~ 30% of 
the level of the Galactic emission at the same energy. 



Table 5. Best-fit parameters using a model with redshifted O lines 



Parameter 


Offsct-A Case 


Offset-B Case 


ICM kT (keV) 


c 41+0.28 


6.02 ±0.22 


ICM Z (solar) 


0.20 ±0.04 


0.21 ±0.04 


ICM z 


0.1756 (fixed) 


0.1756 (fixed) 


ICM Normalization 11 


2.08 ±0.03 x 10 -5 


2.00 ±0.03 x 10 -5 


OVII E (eV) 


488.22 (fixed) 


488.22 (fixed) 


O VII I b c 


(<3.7x 10 -8 ) 


(<3.9x 10 -8 ) 


O Villi E (eV) 


555.99 (fixed) 


555.99 (fixed) 


O VIII I h c 


4.0 (< 13.4) x 10~ 8 


6.8 (< 15.8) x 10~ 8 


FI/BI ratio 


0.89 ±0.01 


0.90 ±0.01 


Ax'/Adof 


0.78/2 


2.25/2 



^ J n c !iH (1V/4:Tt(Da(1 + z)Y per solid angle m units of 

10 cm -5 arcmin -2 , where n c is the electron density, njj the 

hydrogen density, and Z?a the angular size distance. 

b In units of photons cm -2 s -1 arcmin -2 

c Upper limits are quoted in 2a confidence level. 



4. Discussion 

We have observed A2218 with the XIS instrument onboard Suzaku to search for the redshifted O emission lines from 
the WHIM, which is possibly forming a large-scale filament around the cluster. The object was selected because its 
rcdshift (z = 0.1756) would allow the XIS to separate the WHIM lines from the Galactic (z = 0) ones, and because an 
elongated structure in the line-of-sight direction was suggested from the previous studies of this cluster. We detected 
no redshifted O lines, and set a constraint on their intensities. In this section, we compare our constraints on the line 
intensity with those reported in other works and also discuss the future prospect of studying the WHIM emission with 
the Suzaku XIS. 

4-1- Comparison with other results 

Kaastra et al. (2003) and Finoguenov et al. (2003) reported positive detections of O lines around clusters of galaxies 
based on XMM-Newton observations. Kaastra et al. (2003) detected significant O lines in three clusters, Sersic 159-03, 
MKW 3s and A2052, and possible, but not uniquely proven, O lines in the Coma and A1795 clusters. Finoguenov et al. 
(2003) detected the O lines in the outskirts of the Coma cluster, in particular, the Coma-11 field. Fig. 7 compares the 
O line surface brightness I for their observations and our results. The left and right panels are for the O VII and O VIII 
cases, respectively. The intensity of the O lines of Galactic emission measured by a microcalorimeter experiment of 
McCammon et al. (2002), the one compiled by Lumb et al. (2002) and the one estimated from our Offsct-A observation 
are also shown. The surface brightness quoted for Kaastra et al. (2003) was calculated using the temperature and 
emission measure in their Table 7 and the metal abundance from Tables 4 and 5 of Tamura et al. (2004), assuming the 
solar abundance ratio Fe to O (Fe/O number density ratio of 0.55; Anders & Grevesse, 1989). The surface brightness 
levels of the O VII or O VIII emission lines that were reported so far from the XMM-Newton observations are similar to 
or higher than that of the Galactic emission. In contrast, the upper limit of O VII line intensity in the A2218 vicinity 
obtained in this work is about six times lower than the Galactic level. The upper limit for O VIII is also lower than 
the level reported as a positive detection in other works. The tight upper limit we obtained here demonstrates the 
good spectral capability of the XIS, in particular below 0.5 keV. 
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Fig. 6. A2218 BI (black) and FI (red) spectra and models with redshifted O lines. The two emission lines indicated with solid 
lines in the 0.4-0.6 keV range show the maximum allowed level (2<r) of Ovu and Oviii lines. Dashed lines indicate the back- 
ground emission, (upper left) without considering systematic uncertainty, using the best-fit model of Offset-A observation as the 
Galactic emission, (upper right) without considering systematic uncertainty, using Offset-B observation, (bottom) the model of 
maximum O line intensity after considering systematic uncertainties: with ARFs of 20% thicker contaminant and 10% fainter 
Galactic model than the one determined from the Offsct-B observation. The O VII and O VIII surface brightness was 1.1 X 10" 7 and 
2.7 X 10~ 7 photons cm -2 s _1 arcmin -2 , respectively. 



4-2. Constraint on the WHIM density 

If we assume that the O VII emission line is produced in a cloud with uniform density and temperature under 
collisional ionization equilibrium, the surface brightness / is determined by the electron density n e , the hydrogen 
density tch, path length L, and metal abundance Z of the cloud as 

I = C{T) (l + z)- 3 n c n H ZL, (1) 

where C(T) is a coefficient that depends on the temperature of the cloud. Here, we will constrain the density of 
the cloud assuming the temperature to be T ~ 2 x 10 6 K, since the OVII line is strongest at this temperature. 
Substituting C(T = 2 x 10 6 K) with the value calculated with the SPEX code (Kaastra et al., 1996), our constraint, 
/ < 1.1 x 10~ 7 photons enr 2 s _1 arcmin - at z = 0.1756, and electron-to-hydrogen number density ratio n e /nn = 1.2 
for ionized gas, gives the following condition 

-1/2 , j ,-1/2 



0.1 z Qj 

Note that L = 2 Mpc corresponds to 11', the average diameter of the annular region where we extracted the spectra. 
The overdensity 5 = nn/nn of this cloud is 

-1/2 / j s-1/2 

' 2Mp^J ( 3 ) 



0.1 Z, 
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OVII 
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A2052 Lumb A221E 

MKW3s McCammon Offset-A 



Fig. 7. Comparison of Ovn (left) and Ovin (right) surface brightness. From left to right, those in the Coma-11 field (Finogucnov 
ct al., 2003), Sersic 159-03, MKW 3s and A2052 (Kaastra ct al., 2003), Galactic emission of McCammon et al. (2002), Lumb et al. 
(2002), and our Offset-A observation, and the upper limits in A2218 outskirts (this work) are plotted. 

where n-R = Xil\,p cr i t (l + z) 3 /m p = 1.77 x 10 -7 (l + z) 3 cm -3 is the mean hydrogen density in the universe, in which X = 
0.71 is the hydrogen-to-baryon mass ratio, fib = 0.0457 is the baryon density of the universe. p CT n = 9.21 x 10 -30 g cm~ 3 
is the critical density of the universe, and m p is the proton mass. Even though this level of overdensity is much higher 
than the typical WHIM density (S ~ 10), it shows that Suzaku can certainly detect the high-density end of the WHIM 
distribution that is predicted to exist near clusters. 

4-3. Prospect for the WHIM observation with the XIS 

Although we have set a tight constraint on the intensity of redshiftcd O lines, the inherent sensitivity of the XIS was 
presumably not achieved, since the redshiftcd OVIII line fell almost on the Galactic OVII line, while the redshiftcd 
O VII line overlapped with the instrumental neutral O edge. Therefore, better sensitivity is expected for clusters with 
a redshift such that the redshiftcd lines arc free from the Galactic and instrumental features. For example, a redshift 
of z ~ 0.07 or z > 0.2 is suitable for this purpose. The former case is more promising, because the surface brightness of 
high-z clusters is lower and because the contamination of bright ICM emission in the center of the cluster makes the 
detection of warm-hot emission in the outskirts difficult for distant clusters. Note that the XIS is capable of separating 
a redshift difference of Az ~ 0.07 at the O line energy. 

Because of the contamination on the XIS OBF, the effective area at the redshiftcd O VII line had dropped by ~ 25% 
at the time of the A2218 observation. If there had been no contamination on the OBF, we could have obtained more 
photons by this fraction. This contamination also caused a large systematic uncertainty corresponding to about 5% 
of the OVII line flux. Since the thickness and chemical composition of the contaminant will be better understood 
as more data are accumulated, the uncertainty will be less influential in future observations. Furthermore, we hope 
that the original XIS sensitivity in the soft X-ray energy range can be recovered in the future by hardware changes 
including warming the contaminated OBFs. 

It should be noted that in most cases concerning the search for faint WHIM emission with CCD sensors, the variation 
of Galactic emission causes a larger uncertainty than the statistical one. Offset observations, such as those carried out 
in our work, are always desirable to determine the background emission level reliably. Looking at a somewhat longer 
time frame, observations with microcalorimeters are promising, because their outstanding energy resolution (AE < 7 
eV) for extended sources can easily separate the WHIM emission from the Galactic spectral features. For example, two 
orders of magnitude higher sensitivity than the present work is expected with a small mission such as DIOS (Diffuse 
Ionized Oxygen Surveyor; Ohashi ct al., 2006), which is proposed to be launched in early 2010s (Yoshikawa et al., 
2004). 

5. Summary 

We used Suzaku XIS observations to constrain the intensity of O emission lines around the cluster of galaxies 
A2218. After considering systematic uncertainties mainly caused by the Galactic emission and the uncertainty in 
the detector response, we obtained upper limits for the surface brightness of the OVII and OVIII lines of 1.1 x 10~ 7 
photons cm~ 2 s _1 arcmin~ 2 and 3.0 x 10~ 7 photons cm -2 s _1 arcmin -2 , respectively These upper limits are signifi- 
cantly lower than the previously reported fluxes around other clusters of galaxies. Our tight constraints demonstrate 
the sensitivity of the XIS for redshiftcd O lines. 
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